Density and spin response functions in ultracold fermionic atom gases 
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We propose a new method of detecting the onset of superfluidity in a two-component ultracold 
fermionic gas of atoms governed by an attractive short-range interaction. By studying the two-body 
correlation functions we find that a measurement of the momentum distribution of the density and 
spin response functions allows one to access separately the normal and anomalous densities. The 
change in sign at low momentum transfer of the density response function signals the transition 
between a BEC and a BCS regimes, characterized by small and large pairs, respectively. This 
change in sign of the density response function represents an unambiguous signature of the BEC to 
BCS crossover. Also, we predict spin rotational symmetry-breaking in this system. 
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The BEC-to-BCS crossover has drawn renewed inter- 
est in recent years due to experimental progress in atom 
trap systems Q- Several groups j2| achieved a Bose- 
Einstein condensate (BEC) in which the fermions form 
non-overlapping (Shafroth) pairs Q in a two-component 
Fermi gas. Despite considerable effort, much work is still 
needed to understand the other limit in which the pairs 
are large, and the pairing occurs in momentum space, 
similarly to the Bardcen-Coopcr-Schrieffcr (BCS) state 
of superconductivity in normal metals. Questions still 
remain regarding whether or not the system is super- 
fluid when passing through the crossover with tempera- 
tures of the order of twenty percent of the Fermi temper- 
ature. The appropriate characterization of the system 
is still open to debate. Recently, several groups have 
claimed 0, [3 to have reached the supcrfluid state on the 
negative scattering length side of the Fcshbach resonance. 
Despite evidence of fermionic pairing, these claims are 
still subject to intense discussions as no definitive proof 
of superfluidity |^ is available yet. 

Vivcrit et al. have suggested recently that the shape 
of the atomic momentum distribution at low tempera- 
tures is very sensitive to the sign and size of the scat- 
tering length. Altman et al. have proposed to utilize 
density-density correlations in the image of an expanding 
gas cloud to probe complex many-body states of trapped 
ultracold atoms. Also, Bruun and Baym showed that 
scattered light from a fermionic gas would exhibit a large 
maximum below the superfluid critical temperature and 
therefore it can be used to detect the superfluid transi- 
tion. However it is not clear how this behavior changes 
close to the crossover, where the actual experiments are 
performed. 

In this paper we propose a new diagnostic method for 
fermionic atomic gas condensates. By studying the zero- 
temperature evolution of the density and spin response 
functions, as a function of the scattering length of the 
interaction, we show that the density response function 



changes sign across the BEC to BCS crossover and that 
this change can be used to experimentally distinguish the 
BEC from the BCS state of the system. 

Our model system consists of fermionic atoms in two 
hypcrfinc states interacting via a flnite-range attractive 
interaction. Based on a variational approach, we also 
derive a sum-rule satisfied by the spin-spin correlation 
function at g = 0. We find that the density and spin re- 
sponse functions are given by the sum and the difference 
of independent contributions arising from the normal and 
anomalous density, respectively. By measuring the mo- 
mentum distribution of the two response functions it is 
possible to infer separately the normal and anomalous 
densities. Finally, we predict symmetry breaking of the 
spin-rotation invariance. This prediction is characteris- 
tic to our model and, therefore, provides an experimental 
check for its applicability to real systems. 

We begin by considering the two-body Hamiltonian 



H 



^fe,' ^fc^ 



-V -at at 

2 ^q-,im,nj O-km "^f 



i ^P-qj ^k+qn , (1) 



where {aj^^,afei} are the particle creation and annihila- 
tion operators corresponding to single-particle states of 
linear momentum k and fermion type i. The atomic lev- 
els, which we associate with the fermion types, are eigen- 
states of the Hamiltonian of an ion (with integer nuclear 
spin /) interacting with an electron (spin s = ^) 
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where A denotes the strength of the hyperfine interac- 
tion and B is the magnetic field, while fi^ and denote 
the electron and nuclear magnetic moments, respectively. 
We introduce the total angular momentum, F ~ I + s. 
The total angular momentum projection, Mp, is the only 
good quantum number at finite B. The atomic spectra 
considered in our model are depicted in Fig. 

In the Hartree-Fock-Bogoliubov (HFB) formalism, one 
first introduces the quasi-particle creation and annihila- 
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FIG. 1; Atomic states involved in the Feschbach resonance at 
small mae netic field, B, for (a) ®Li and (b) *°K. The hyperfine 
couplings between states allowed by the selection rules are 
represented by dotted lines. In ^Li, the bound state is formed 

"^^K, it is formed from 



from the -F = 3/2 states while, 
1 1 — I) and the lowest eigenstate 
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tion operators, {(3ki, in terms of the particle creation 
and annihilation operators, via the Bogohubov-Valatin 
transformation 0. Then, the ground state, |<i>), of the 
Hamiltonian JQ), is obtained via a variational ansatz, 
such that /?fci|$) = 0, for all {ki} labels. We will refer to 
this state as the HFB wave function. 

While the above approach is general for any arbitrary 
multi- level fermionic Hamiltonian ^3 j for the purpose of 
the present discussion we will confine ourselves to a two 
level (one-channel) model, such as discussed in [ill IT^. 
In this model the HFB wave function has the BCS form 
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subject to the normalization condition. 
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For zero magnetic field, B, the states | t) and | J,) refer 
to the two fermion types \F Mp) = ||,|) and ||,— 5) 



we have 



for ^Li, and ||, — |) and ||, — |) for "'"K. For finite values 
of B 
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TABLE I: Parameters entering the modified Pauli matrices. 
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for °Li, while for we have 
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The parameters (a, 6, c, d) are related to the hyperfine 
mixing angles, such that 

a = sin 01,6 = cos 01 , c = sin 02 , (i = cos 02 , (7) 

for ^Li, and 

a = cos 01,6 = sin 0i , c = cos 02 , d = sin 02 , (8) 

for ^"K. When the magnetic field is zero, the above pa- 
rameters are given by the appropriate Clebsch-Gordan 
coefficients, while for large fields we retrieve the unmixed 
phase, as illustrated in Tabled 

The ground-state properties are described by the nor- 
mal and anomalous densities defined as 

Pk = (<i>|a|,^akTl<i>> = bkP, (9) 
Kk = (^la.kiakTl^") = t'kUk, (10) 

while the mean total particle-density of the system, po, 
is given by 
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The ground-state ansatz Q provides a smooth inter- 
polation between the BCS and the BEC regimes. We 
have recently used this model to study the prop- 
erties of the BEC-BCS crossover, for a short- (but fi- 
nite) range and attractive interaction. In the dilute limit 
the model is equivalent to the zero-range (contact) inter- 
action Hamiltonian, initially discussed by Leggett [l3| . 
The study presented in Ref. was carried out by 
modifying the scattering length of the interaction, while 
keeping the density and range of the interaction fixed. 
Figure 12 illustrates the momentum distributions of the 
normal and anomalous densities. The mean-field solu- 
tion predicts the crossover occurs when the minimum 
in the quasi-particle energy spectrum shifts from a fi- 
nite (BCS) to zero-momentum value (BEC). As such, 
the presence/absence of the singularity in the momentum 
distribution of the density of states represents a unam- 
biguous signature of the crossover (see inset in Fig.l). 

Here, we focus on the changes in the two-point corre- 
lation functions as we evolve the system from the BEC 
to the BCS regimes. We begin by defining the response 
functions 
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FIG. 2: (Color online) Momentum distribution of the nor- 
mal and anomalous densities, as a function of the scatter- 
ing length, at fixed system density, po (kpir) ~ 0.37) . 
Inset shows the density of states. We use the notation 
1] = {aokpY^- 



where p[t) and Sfi{v) denote the particle- and spin- 
density operators 

p(r)=^al(z|<5(r-ri)|j)a,, (14) 
5,,(r) = i^al(*|a^<5(r-ri)U)a,. (15) 

Here, denote the Pauli matrices. Since the single- 
particle states are plane waves, i.e. \i) = e^^^ '^Xi , where 
X denotes the (| t), | J,)) spinors, then we can calculate 
the particle-density matrix element as 



(i|<5(r-ri)|j) = 5.,, eHk.-k.)-r ^ 
and the spin-density matrix element as 
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The density response function in (|12|l is a scalar, while 
the spin density response function in (|13|l is a tensor. 
For the one-channel model the spin-response tensor is 
diagonal, and, provided that the Hamiltonian and the 
total spin operator commute, the diagonal components 
are equal. Then, the density and spin response functions 
correspond to the sum and difference of separate normal 
and anomalous density contributions, respectively. The 
nonlocal response functions can be written as 



Sp[q) = ^Kiq) - lp{q) , 
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FIG. 3: (Color online) The density and spin instantaneous 
response function for an electron gas, as a function of the 
scattering length (or r) — (aokp)'^), at fixed density, po. 



Equations H18() and H19|l represent a general result for 
the HFB mean-field approximation of the ground-state 
of the multi-level Hamiltonian For the one-channel 
model, and disregarding for now the hyperfine nature of 
the atomic levels involved (similarly to the case of an 
electron gas), the normal and anomalous density contri- 
butions to the response functions arc: 
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The response functions Sp{q) and Sa{q) arc shown in 
Fig. 13 The density response function, Sp{q), changes 
sign at low momentum, as the bound state disap- 
pears, while the spin response function, S^{q), changes 
smoothly from the BEC to the BCS regime. In the BCS 
limit, the condensate wave function, Kk, appears in a very 
narrow region around the Fermi momentum (see Fig.lJJ. 
Since the width of this distribution tends to zero in the 
limit ao — > —0, so does the anomalous density contri- 
bution, Xk((7). In turn, the density and spin response 
functions will be equal in this limit. 

The mean-field approach shows that a measurement of 
the density response function is a signature of the BEC- 
BCS crossover. In the dilute limit, the crossover coincides 
with the singularity in the scattering length, which in 
turn corresponds to the change in sign of the density re- 
sponse function. Figure2](top) shows the dependence of 
Iniq) and 1p{q) for g = 0, as a function of the scattering 
length ao at fixed density po (kpir) ~ 0.37). The corre- 
sponding density response function Sp{0) is depicted in 
Fig. ^ (bottom) . The density response function changes 
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FIG. 4: (Color online) Zero-momentum transfer normal and 
anomalous density contributions, Ik{0) and lp{0), together 
with the density response function Sp{q) for q = 0. 



of the spin response function is broken. This symmetry- 
breaking effect is a prediction of our model, and arises 
as a consequence of the fact that our effective interaction 
involves a restricted set of atomic levels. 

In conclusion, in this paper, we have shown that much 
information about the crossover regime can be gained by 
experimentally studying the density and spin response 
functions. Within the framework of the mean-field re- 
sults, we show that the normal and anomalous densities 
can be accessed from the momentum distribution of the 
response functions. The spin response function changes 
smoothly across the crossover, while the density response 
function changes sign, and thus represents a signature 
of the crossover. The spin response at zero-momentum 
transfer satisfies a sum rule, and is sensitive to the in- 
teraction. When taking into account the hyperfine struc- 
ture of the interacting levels, our model predicts that the 
rotational invariance of the spin response, normally asso- 
ciated with an electron gas subject to a spin-independent 
interaction, is broken. 



sign between the BCS and the BEC limits. 

The one-channel model predicts that the spin response 
function at zero momentum transfer, S^iq = 0), is in fact 
a sum rule, i.e. 
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or Tp{0) ^-^^(O) = 1 is independent of ij ~ (aokp)"^, as 
shown in Fig. 0] (top). The above is obtained by using the 
definitions © and (|10|l , together with the normalization 
condition, Eq. Q). 

We now consider the modification of the electron gas 
results, due to the hyperfine nature of the interacting 
atomic levels. We find the spin response function is pro- 
portional to the "weighted" difference of the normal and 
anomalous density contributions derived earlier, Ik (9) 
and 1p{q)- The weighting factors depend on the atom 
specie in the system, and the atomic levels involved in the 
interaction. Since the spin response function is only sen- 
sitive to the electron spin operator, then in order to find 
the modification of the spin response function, we need 
to calculate the matrix elements (xilcr^lxj) in Eq. (|17|l . 
The parameters entering the modified Pauli matrices are 
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and limiting values as a function of the magnetic field, 
B, are listed in Table We note that in the large field 
limit, the spin response in a ^Li fermionic atom gas has 
opposite sign as compared to the case of a fermionic 
atom gas. Irrespective of the actual B value, we obtain 
that for a fermionic atom gas, the rotational symmetry 
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